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SUMMARY 


The  first  step  in  the  development  of  a hybrid  com- 
puter simulation  of  the  data-based,  5 DOF,  3K-SES  program 
that  will  operate  in  real-time  has  been  completed.  A 
3 DOF  (Flat)  turn  model  for  an  approximate  and  simplified 
XR-3  craft  was  used  to  provide  the  equations  of  motion  in 
order  to  become  familiar  with  1)  the  turning  character- 
istics of  the  CABSES,  2)  the  requirements  of  the  (NPS)  NAV 
PG  School  hybrid  computer  with  regard  to  signal  conversion, 
timing,  accuracy  and  loop  delay  due  to  both  signal  con- 
version and  (right-hand-side)  computation  and  3)  the 
approximate  size  of  program  that  can  be  computed  on  the 
existing  hardware  of  the  NPS  facility. 

In  parallel  with  this  investigation,  the  5 DOF  data- 
based  digital  program  for  the  3K-SES  has  been  converted 
from  CDC  to  IBM  FORTRAN  and  is  now  operational  and  will 
serve  as  the  standard  for  evaluation  of  the  final  hybrid 
simulation  for  the  3K-SES . 

Results  of  the  hybrid  computer  simulation  of  the  3 DOF 
turn  model  of  the  XR-3  indicate  that  the  installation  of  the 
5 DOF  program  of  the  3K-SES  should  proceed  (in  the  next 
phase)  as  a straight  forward  application  of  the  right-hand- 
side  computations  within  the  digital  computer  and  addi- 
tional degrees  of  freedom  on  the  analog  computer.  The 
successful  operation  of  the  hybrid  computer  simulation  will 


4 


depend  upon  ability  of  the  10  year-old  hardware — in 
particular  the  digital  computer — to  compute  the  many 
computational  operations  required  for  the  right-hand- 
side  of  the  equations  of  motion  in  the  time  frame  necessary 
to  approximate  real-time  manual  control. 
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Reference  1 presents  the  Rohr  Marine,  Inc.  analysis 
of  the  stability  and  maneuverability  of  the  3K-SES  under 
direct  manual  control.  The  recommendations  contained  in 
Reference  1 included  a suggestion  to  continue  the  investi- 
gation of  manual  control  using  a more  realistic  model  of 
the  response  of  the  human  operator  to  emergency  conditions. 

If  a real  time  simulation  can  be  developed  it  can 
provide  a means  of  inserting  an  actual  human  operator  in 
the  loop,  as  well  as  any  selected  model  of  the  operator. 

A hybrid  realization  of  such  a real  time  simulator  appears 
to  be  feasible  and  would  provide  a convenient  hands-on 
facility  for  a variety  of  human  operator  studies,  plus 
many  other  studies  in  control  and  guidance. 

This  report  presents  the  results  of  the  first  phase 
of  a feasibility  study,  in  which  it  is  planned  to  develop 
a real  time  simulation  using  the  NPS  Hybrid  Computer 
which  consists  of  a CI-5000  analog  unit  and  an  XDS-9300 
digital  unit.  The  initial  studies  reported  here  are 
concerned  with  development  of  3 DOF  models  for  the  XR-3 
and  the  3K,  as  a preliminary  to  the  development  of  a 5 DOF 
Hybrid  model. 


- 
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IIA.  THE  3 DEGREE  OF  FREEDOM  MODEL 

1.  INTRODUCTION 

The  development  of  any  computer  model  requires  famili- 
arity with  both  the  system  equations  and  the  computer  to  be 
used.  In  addition,  one  needs  quantitative  data  to  insert 
the  model  into  the  computer,  and  one  also  needs  typical 
performance  data  such  as  experimental  test  results  in  order 
to  validate  and  verify  the  model.  To  develop  a hybrid  model 
of  the  3K-SES,  we  chose  to  start  by  modelling  the  ”R-3, 
because  we  have  the  required  quantitative  data,  we  have  a 
good  6 DOF  digital  program,  and  we  have  good  experimental 
data,  plus  availablility  of  the  XR-3  craft  for  test  as  needed. 
We  also  decided  to  start  with  a 3 DOF  model,  for  several 
reasons : 

a)  It  is  the  simplest  model  which  will  permit  turning 
tests  that  are  verifiable. 

b)  The  complexity  of  the  initial  model  will  thus  be 
minimal,  but  we  can  readily  add  two  more  degrees 
of  freedom  with  confidence  in  the  resulting  5 DOF 
model,  because  the  similarity  of  the  differential 
equations  will  permit  ready  application  of  the 
techniques  we  develop  in  perfecting  the  3 DOF  model. 

In  addition  to  this  decision  to  start  with  a 3 DOF  model, 
we  chose  to  simplify  the  model  by  simplifying  the  right-hand- 
side  of  the  equations.  Since  the  terms  on  the  right-hand- 
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aide  are  arithmetic  calculations  to  be  done  in  the  digital 
unit,  and  because  the  many  terms  are  additive , truncation 
of  the  right -hand- side , with  simplification  of  the  terms 
retained,  should  have  little  effect  on  results  except  some 
decrease  in  quantitative  accuracy.  This  permits  us  to  work 
with  a less  complex  model  in  the  initial  stages,  and  addition 
of  terms  at  a later  point  should  be  a simple  matter.  It 
will  also  have  the  advantage  of  permitting  us  to  assess  the 
effects  of  adding  terms  on  the  computational  delays  and  on 
the  limits  of  capability  of  our  hybrid  installation. 

2.  EQUATIONS  OF  MOTION  OF  THE  3 DOF  MODEL 

The  three  degrees  of  freedom  to  be  used  in  this  devel- 
opment are:  SURGE,  SWAY,  YAW.  As  a result  the  following 
assumptions  and  approximations  are  needed: 

1.  Roll  is  2ero 

2.  Pitch  and  draft  remain  constant 

3.  The  origin  of  the  ship’s  coordinate  system  is 
located  at  the  center  of  gravity. 

4.  At  t = 0 the  ship  coordinate  system  and  the  fixed 
(geographical)  coordinate  system  have  the  same 
origin  and  the  same  coordinate  axes. 

5.  At  t * 0 the  ship  is  heading  in  the  positive 

X_ ...  direction, 
nav 


11 


/'9ur*  '■  *»"*»  op  coo**,* 


SVSTEMS. 


ship's  path 


„ ORIGIN  OF 

SH,P#  coordinates 


* V 

K*  v * 


\ I 

\ -y  | 


£ _ u - 


ORIGIN 
OF 
. F;X£D 

ars1 


-a  f/  1 ' 

3 i \ 


/ i 
/ i 

i 

i 

i 

i 

i 


aSSa- 


OrGhr  rw)N 


12 


The  equations  of  motion  can  then  be  approximated  by: 


SURGE 

m( u - 

vr)  ® 

X 

(1) 

SWAY 

m(v  + 

ur)  * 

Y 

(2) 

YAW 

• 

rz  r 

m 

N 

(3) 

Navigation 

xQ  ■ u cos  Y - 

v sin 

Y 

(4) 

31  u sin  Y + 

v cos 

Y 

(5) 

v3  .-f? 

♦ v2 

(6) 

S ■ tan 

-1  (Iv_ 
u 

•) 

(7) 

In  the  above 

equations : 

m « mass 

of  the 

rigid 

ship 

I * moment  of  inertia  about  the  z-axis 
2 


u 


vg  cos  8 = velocity  in  the  x-direction  (SURGE) 


v » V sin  8 * velocity  in  the  y-direction  (SWAY) 
s 

r * f * angular  velocity  about  the  z-axis 
Y * heading  (yaw)  angle  (see  Fig.  1) 

3 * drift  angle  (see  Fig.  1) 

>5  » rudder  angle,  + 5 produces  a left  turn  (see  Fig.  1) 


Also: 


:o  mJ 

fkQdt 

u = J\ adt 

r -/rdt 

= j 

(y  dt 

v 3 /v dt 

Y = / Ydt 

O J 

•'0 

J 

J 
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The  right-hand-sides  ( RHS ) of  the  equations  of  motion 
have  been  designated  X,  Y,  N.  These  symbols  represent 
summations  of  forces  (X,Y)  and  moments  (N) , and  each  summation 
may  consist  of  many  terms.  For  the  initial  simplified  model 
we  have  discarded  all  terms  in  the  summations  except  the 
few  that  are  obviously  major  contributions.  The  terms  re- 
tained are  considered  sufficient  to  produce  a working 
simulation  which  should  provide  qualitative  agreement  with 
reference  data.  At  a later  point  the  model  may  be  refined 
by  adding  terms  and  adjusting  terms  and  adjusting  coefficients. 

After  reduction  of  the  right-hand-sides,  these  force 
summations  become: 

X - GT  cos  i - CDXu2  (8) 

Y - GT  sin  5 - C^vM  (9) 

N - -(GT  sin  mQ  + vlvllu  (1Q) 

where:  (see  Fig.  2 also) 

GT  is  the  gross  thrust  in  pounds 

5*5  =5  * angle  of  all  thrust  vectors 

3 p 

CQX  ■ coefficient  of  drag  in  the  direction  of  surge 
* coefficient  of  drag  in  the  direction  sway 

1 (see  Fig.  2)*  thrust  lever  arm 
o 

1^  (see  Fig.  2)*  sway-drag  lever  arm 

In  further  explanation  of  equations  8,  9,  10,  note 
that  we  have  chosen  to  model  the  drag  forces  as  lumped 
parameters  acting  at  two  points  P and  p'  (see  Fig.  2) . 


1 


n 

i 


u 


STARBOARD 


NOTE  : 8 «S  -8 
* P 

P - P1  ARE  EQUIVALENT  FORCE  CENTROIDS 
Figure  2.  FORCES  ACTING  FOR  A RIGHT  TURN. 


The  length  1 is  the  distance  from  the  stern  to  the  center 
of  gravity,  and  1^  is  the  component  distance  (in  the  x- 
direction)  from  the  center  of  gravity  to  the  points  P-P'. 

1Q  is  a fixed  parameter  for  our  initial  studies  (it  could  be 
changed  by  ship  loading)  but  1 is  an  adjustable  distance 
which  we  will  use  as  a parameter  to  make  the  model  fit  test 
data. 

The  coefficient  CQX  is  a lumped  parameter  equivalent 

for  the  skin  friction  effect  (slender  body  theory) . The  term 
2 

CDXu  models  the  hull  drag,  aerodynamic  effects  and  hydro- 
dynamic  effects.  The  bubble  drag  is  not  included  in  the 
initial  model  (will  be  addded  later)  because  our  initial 
efforts  are  aimed  at  evaluation  of  our  hybrid  structure, 
and  a simple  model  above  critical  speed  will  suffice  for 
this.  In  like  manner,  in  eqn.  9,  the  coefficient  CQy  models 
the  cross  flow  drag  effects.  Note  that  the  added  mass  effect 
has  been  temporarily  omitted  since  its  effect  is  also  of 
second  order. 

Substituting  equations  8,  9,  10  into  equations  1,  2,  3 
provides  the  basic  equations  of  motion  for  the  3 degree  of 
freedom  model: 

u =»  (GT  cos  6 - Cnv u2)  + vr  (11) 

m Da 

v =»  -i-  (GT  sin  5 - Cj^v|v|)  - ur  (12) 

r * -p-  ^-(GT  sin  <S)  1q  + (C^v | v|  ) 1^  j (13) 

Equations  11,  12,  13  are  applicable  to  any  SES  craft. 
Selection  of  parameter  values  to  be  used  in  the  equations 


16 


f 


would  depend  on  the  particular  craft,  but  the  procedures 
used  to  identify  the  lumped  parameter  values  are  independent 
of  the  craft  type. 

3.  IDENTIFICATION  OF  CRAFT  PARAMETERS 

In  order  to  evaluate  the  parameters  needed  for  the  equa- 
tions, the  craft  must  be  available  for  measurement and  test 

(or  suitable  data  must  be  available) . The  parameters  m,  I , 

z 

1Q  are  constants  which  must  be  known,  and  GT,  5,  u,  v,  r 
can  be  measured  in  steady  state  (if  needed) . Then  C„v,  1 , 

DA  u> 

C^,  can  be  calculated  as  is  shown  below: 

Dy 

A.  Surge  Drag  Coefficient,  CDX 

The  craft  is  run  on  a straight  course  with  GT  = 
constant  and  given  rudder  conditon  5 . When  steady  state  is 
reached  u is  measured.  Then 

u . 0 - -i-  (CT  cos  « - CpjjU J3)  (14) 

from  which 

CDX  = GT  cos  5/u gS  (15) 

where 

GT  = gross  thrust  in  pounds 

u = steady  state  speed  in  feet/sec 
s s 

5 = rudder  angle 

B.  Sway  Drag  Moment  Arm,  ly 

Let  GT  = constant  and  5 * 6^,  a fixed  angle.  When 
the  craft  reaches  steady  state  in  the  turn  the  sway  equation 
becomes 

v = 0 = (GT  sin  5 - C^vj  v | ) — ur  (16) 
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from  which 


C_  v = GT  sin  { - m u r 
Dy 

The  steady  state  yaw  moment  equation  becomes 


1 

L“ 


(GT  sin  5)lo  + (CpyVl v| ) lo 


from  which 


CDyv2  = (GT  sin  5) (lo/l  } 

^ i(ll 

Combining  equations  17  and  19: 


GT  sin  6(1  - *o/  ) 

0) 


or  if  we  wish 

1 9? -s  -n 5 (21) 

W xo  GT  sin  6 - m u r 

From  Fig.  2 and  eqn.  20  it  is  clear  that  the  turning 

rate,  r,  is  a function  of  the  distance  1 , with  all  other 

01 

parameter  values  fixed.  Since  it  is  customary  to  command 

a turning  rate,  if  we  choose  r we  can  then  evaluate  1 . For 

01 

the  XR-3,  using  steady  state  data  obtained  with  our  6 DOF 
simulation, 

r = -0.0266(1  - 12/j_  ) radians/sec  (22) 

0) 

A plot  of  this  relationsip  is  shown  on  Fig.  3. 

C.  Sway  Drag  Coefficient,  C_ 


Again  assuming  steady  state  conditions  with  fixed 

thrust  and  fixed  rudder,  the  sway  equation,  eqn.  17,  provides: 
_ _ -m  u r + GT  sin  6 
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Figure  4-B  FLAT  TURN  DYNAMICS  OF  THE  XR-3,  PLOT 
OF  X ve.  Y. 

r*  .33  RAD/SEC  (fl  * I ft) 


Figure  5-B  FLAT  TURN  DYNAMICS  OF  THE  XR-3,  PLOT 
OF  X vs.  Y. 

r = O.l  RAD/SEC  ( 0 = 2.5  ft) 
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For  a given  turning  rate,  r,  we  know  all  values  except  v, 
which  can  be  calculated  using  the  surge  equation,  eqn.  11 


to  provide: 


CQXu  - GT  cos  5 
m r 


(24) 


4.  DIGITAL  SIMULATION  OF  THE  SIMPLIFIED  XR-3  EQUATIONS  IN 
3 DEGREES  OF  FREEDOM 

The  equations  of  motions  as  previously  derived  were 

programmed  into  the  IBM  360/67  computer.  A fixed  step  size 

was  used  to  verify  the  program  and  validate  the  model. 

Using  a turn  rate,  r = .33  rad/sec,  and  an  integration  step 

of  0.1  seconds  the  simulation  was  run  for  120  seconds  of 

problem  time,  then  repeated  for  an  integration  step  of  0.01 

seconds.  Comparison  of  results  showed  that  an  integration 

step  of  0.1  seconds  was  suitable,  all  variables  in  the  two 

runs  agreeing  to  within  1% . 

Figure  4A  and  4B  show  plots  of  u vs  time  and  x vs  y 

for  a flat  turn  with  r = .33  rad/sec  (1  = 1 ft.).  Figures 

5A  and  5B  show  similar  plots  for  r = 0.1  rad/sec  (1  = 2.5  ft.). 

0) 

A number  of  additional  runs  were  made  but  results  are  not 
given  here.  Based  on  these  runs  we  feel  that  the  program 
has  been  satisfactorily  verified,  i.e.,  it  is  solving  the 
equations  correctly. 

The  following  observations  seem  appropriate: 

1.  Eqn.  22  and  Fig.  3 are  at  least  qualitatively 

correct. 

2.  The  effect  of  added  mass  is  neglegible. 


i 


r 


1 
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3.  We  have  gained  confidence  that,  given  the  meas- 
ured perfomance  of  any  SES  in  an  approximate  flat  turn,  and 
given  values  for  the  parameters  GT,  5,  u,  r,  m and  1Q,  we 
can  determine  a value  of  1 which  provides  a simulation 
matching  the  measured  performance  reasonably  closely. 

4.  If  the  simulation  is  performed  with  an  initial 
step  in  the  thrust  angle,  5,  approximate  values  for  the  system 
time  constants  can  be  evaluated  from  the  resulting  data. 

For  the  XR-3  operating  in  a flat  turn  the  time  constants  are: 

Variable  Time  Constant 

u 11.0  sec 

v 0.8  sec 

r 1.2  sec 

6 0.9  sec 

To  validate  the  3 DOF  model,  additional  runs  were  made 
using  the  digital  6 DOFOceanics  program  as  modified  by  our 
NPS  group  for  the  XR-3.  Also  the  3 DOF  model  was  run  with 
3K-SES  data,  and  the  5 DOF  model  of  the  3K-SES  was  run.  For 
the  XR-3  tests  a step  of  rudder  of  0.5  radians  was  suddenly 
applied  with  craft  speed  at  20  knots.  For  the  3K-SES  tests 
and  11°  rudder  was  applied  by  ramping  up  for  1.0  second  and 
holding  at  the  11°  value.  Craft  speed  was  60  knots.  Com- 
arison  of  the  3 DOF,  6 DOF  and  3 DOF  Hybrid  runs  for  the 
XR-3  is  shown  on  Fig.  6.  Correlation  is  not  as  bad  as 
appears  by  inspection.  Turning  diameter  and  time  for  a 360° 
turn  are  almost  indentical.  We  presently  attribute  the 
differences  to  the  fact  that  the  3 DOF  model  does  not  include 
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a term  for  bubble  drag. 

Figure  7 compares  the  3 DOF  run  with  the  5 DOF  run 
for  the  3K-SES.  The  results  are  almost  identical,  despite 
simplifications  in  the  3 DOF.  This  gives  confidence  that 
our  modelling  technique  will  work  well  when  we  hybridize. 


5.  HYBRID  SIMULATION  OF  THE  3 DEGREE  OF  FREEDOM  MODEL  OF 


THE  XR-3 

In  developing  a hybrid  model  various  arrangements  and 


structures  (programming)  are  possible.  The  philosophy 
chosen  for  our  initial  hybrid  model  is  shown  in  the  block 
flow  diagram  of  Fig.  3.  Details  of  the  programming  are 
given  in  Appendix  A for  the  simplified  XR-3  model. 


The  simplified  XR-3  model  was  installed  in  the  hybrid 


tained  by  digital  simulation  showed  almost  identical  per- 
formance, differences  being  less  than  2%. 

In  developing  the  hybrid  simulation  a major  concern 
has  been  the  amount  of  time  delay  introduced  into  the  loop 
by  the  A/D  and  D/A  conversions  and  by  the  time  required  to 
perform  the  digital  computations.  We  have  monitored  these 


delays  in  the  3 DOF  model  with  results  as  follows 


1.  The  fundamental  loop  delay  using  D/A  for  u,  v,  r 


2.  A lower  language  subroutine  exists,  called  ADDA 


which  can  be  used  to  replace  the  Fortran  subroutines 


for  A/D  and  D/A.  This  can  reduce  the  delay  by  25% 


3.  Additional  delays  are  introduced  by  using  A/D 


for  thrust  (GT)  and  thrust  angle  (6)  [76%  increase] 

and  by  use  of  the  Graphics  Signal  (Gated)  conversion 
of  the  variables  [67%  increase]  . Total  loop  trans- 
port delay  using  D/A,  both  A/D  units  and  the  Graphics 
signal  conversion  was  approximately  78  ms. 

We  know  how  to  reduce  the  loop  delay  substantially,  but 
in  view  of  the  system  response  time  this  may  not  be  necessary. 

One  possible  source  of  difficulty  in  the  hybrid  simulation 
is  the  cross  coupling  between  the  equations.  When  operating 
near  steady  state  conditions  during  small  yaw  angle  maneuvers 
the  linearized  equations  indicate  that  the  roots  of  the 
characteristic  equation  may  be  very  lightly  damped. 

Consider  straight-ahead  motion  (thrust)  at  constant 
speed.  The  linear  model  sway  and  yaw  equations  are: 


where 


x2  = 

a22X2 

Sa- 

a32X2 

x'  A 

£&u,  AV 

a22  = 

■ fx(x2) 

a32  S 

1 f2(*2) 

joa3 


Combining  these  equations,  the  characteristic  equation 
of  the  simplified  3 DOF  model  is 

s2  + a22s  + uoa32  ■ 0 

The  locations  of  the  roots  of  this  equation  depend  on  the 
values  of  a22  and  a22'  ^oth  which  are  functions  of  the 
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sway  velocity.  For  small  values  of  sway  velocity  the 
roots  are  complex  conjugates  and  approach  the  origin  of 
the  s-plane  as  the  sway  velocity  approaches  zero.  For 
larger  values  of  sway  velocity  such  as  those  encountered 
in  a turn  maneuver  the  roots  are  real  and  negative,  pro- 
viding a well  damped  response  for  yaw  and  sway  motions. 

This  appears  to  be  similar  to  the  Dutch  Roll  effect  in 
aircraft . 

For  straight  ahead  runs  using  the  digital  program  for 
the  3 DOF  model  and  setting  the  initial  value  of  v to  zero, 
no  sway-yaw  motions  are  observed.  However,  if  we  set  a 
small  initial  condition  on  v a transient  oscillation  is 
observed  in  both  sway  and  yaw  as  shown  on  Fig.  9.  This 
verifies  the  existence  of  the  phenomenon,  and  we  expect 
to  encounter  it  in  the  hybrid  since  the  equations  are  the 
same.  However,  in  the  hybrid  there  is  a possibility  that 
the  roots  may  move  into  the  right  half  s-plane  due  to  either: 

1.  DC  unbalance  in  the  analog  computer  integrators,  or, 

2.  Excessive  loop  delay  from  A/D,  D/A  and  digital 
computations . 

The  first  of  these  possible  causes  can  be  avoided  by 
careful  balancing  of  the  operational  amplifiers  before 
running.  The  effect  of  the  loop  delays  can  be  controlled 
somewhat  by  minimizing  the  number  of  terms  requiring  digital 
computation. 

€.  STATUS  OF  THE  3 DOF  MODEL  STUDIES 
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qualitative  and  quantitative  agreement  have  been  achieved 
with  the  5 DOF  digital  simulation  as  a reference.  Future 
studies  will  be  aimed  primarily  at  the  5 DOF  model  of  the 
3K-SES  except  when  we  wish  to  correlate  computer  data  with 
specially  designed  experiments. 
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I IB.  THE  5 DEGREE  OF  FREEDOM  MODEL 


1.  INTRODUCTION 

Having  developed  a simplified  3 DOF  model  of  the  XR-3, 
we  now  proceed  to  use  the  knowledge  thus  obtained  and  apply 
it  to  the  3K-SES . We  do  not  plan  to  develop  a 5 DOF 
hybrid  model  of  the  XR-3.  Our  procedure  is  to  install 
the  NSRDC  5 DOF  program  for  the  3K-SES  in  our  IBM  360/67 
computer,  generate  reference  curves  with  it,  obtain  any 
numbers  we  need  for  our  5 DOF  hybrid  simulation,  then 
verify  and  validate  the  hybrid  model  of  the  3K-SES. 

2.  THE  DIGITAL  SIMULATION  OF  THE  3K-SES  IN  5 DEGREES  OF 

FREEDOM 

This  program  was  obtained  from  NSRDC.  It  has  been 
written  for  a CDC  machine  and  had  to  be  converted  to  the 
IBM  version  of  Fortran  for  our  computer.  This  has  been 
completed,  the  program  is  operational  and  as  far  as  we 
know  has  no  obvious  "bugs".  We  are  familiarising  ourselves 
with  the  program,  and  exercising  it  to  obtain  data  we  will 
need  in  developing  the  5 DOF  hybrid  model,  and  for  validating 
it. 

Details  of  the  conversion-problems  encountered,  their 
solutions,  etc.  are  discussed  in  Appendix  3. 

3.  HYBRID  SIMULATION  OF  THE  3K-SES 

This  is  under  development  at  the  present  time.  Hope- 
fully we  need  only  insert  a new  set  of  constants  in  the 
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equations  of  the  hybrid  XR-3  model,  but  undoubtedly  other 
adjustments  and  modifications  will  be  needed,  particularily 


in  the  trade-off  of  digital  vs  analog  computation  in  order 
to  perform  a real-time  simulation. 

4 . STATUS 

Progress  has  been  satisfactory.  It  is  too  early  to 
predict  future  developments. 


r 

1 ** 


The  NPS  Hybrid  Computer  consisting  of  the  XDS-9300 
digital  and  the  Cl  5000  analog  computers  has  been  pro- 
grammed to  simulate  the  3 DOF  Flat  Turn  maneuvering 
dynamics  of  the  CAB  type  Surface  Effect  Ship. 

The  XR-3  has  been  programmed  with  simplified  modelling 
of  the  drag  forces  and  thrust-input  force  on  the  hybrid 
computer  for  simple  flat  turn  maneuvers.  An  evaluation 
of  the  hybrid  simulation  was  carried  out  by  comparison  of 
the  time  histories  of  the  velocities  and  the  navigation 
variables  against  those  values  obtained  from  the  3 DOF 
digital  simulation  for  the  same  maneuvers.  In  addition, 
the  XR-3  hybrid  3 DOF  results  were  compared  to  the  6 DOF 
loads  and  motion  output . 

The  result  of  this  comparison  was  that  the  hybrid 
computer  simulation  produced  a time  response  to  a step 
rudder  input  that  produced  a navigational  turning  plot  that 
was  in-between  the  all  digital  3 DOF  and  sill  digital  6 DOF 
trajectories. 

Additionally  the  3 DOF  simplified  RHS  (right-hand-side) 
digital  program  was  operated  for  the  same  turning  conditions 
that  was  used  for  the  5 DOF,  DB5IM5D  Program  and  the 
comparison  indicates  close  agreement  in  the  navigational 
trajectories  for  an  11  degree  vectored  thrust  angle  full- 
turn  maneuver. 
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This  feasibility  study  has  shown  that  the  hybrid 
computer  implementation  for  the  5 DOF  3K-SES  in  real- 
time is  possible  with  the  available  equipment  at  the  NPS 
Computer  Laboratory.  Real-time  computation  is  of  course 
only  approximated  when  using  the  hybrid  computer.  The 
loop  transport  delay  inherent  in  this  simulation  is  due 
to  the  signal  conversion  (A  to  D and  D to  A)  and  the 
necessary  digital  computation  for  the  acceleration  forces 
and  moments.  The  magnitude  of  this  delay  time  has  a fixed 
minimum  value  as  a result  of  signal  conversion  and  a maximum 
value  due  to  righ-hand-side  calculation  required  for  the 
acceleration  added  to  the  fixed  time  due  to  the  signal  con- 
version. Real-time  computation  accuracy  is  determined 
principally  by  the  magnitude  of  this  maximum  value  for  the 
loop  delay. 

The  3K-SES  was  approximated  with  the  simplified  3 DOF 
Flat  Turn  equations  and  the  results  were  compared.  It 
was  shown  that  the  full  maneuvering  trajectories  for  an 
11  degree  thrust  angle  compared  very  closely  with  the  5 DOF, 
DBSIM5D  Program  results.  These  results  indicate  that  only 
a few  terms  on  the  RHS  may  be  all  that  are  required  to  simu- 
late the  3K-SES  for  the  type  of  maneuvering  control  studies 
to  be  conducted  in  the  future  on  the  hybrid  computer. 

The  conclusion  of  this  feasibility  study  is  that  suc- 
cessful implementation  of  a hybrid  computer,  real-time 
simulation  of  the  DBS IMS D Program  has  a very  high  prob- 
ability. The  degree  of  complexity  necessary  for  the  RHS 
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APPENDIX  A.  HYBRID  COMPUTER  PROGRAM  FOR  THE 
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3 DOF  (FLAT)  TURN  EQUATIONS 

A general  schematic  diagram  of  the  hybrid  structure 
used  is  shown  in  Fig.  A-l.  The  listing  of  the  digital 
program  for  the  XDS-9300  computer  is  also  included  in  this 
appendix,  and  the  analog  computer  connections  are  shown 
on  Fig.  A-2. 

As  may  be  seen  from  the  listing  of  the  digital  pro- 
gram, the  right-hand-sides  (RHS)  of  the  acceleration 
equations  for  the  3 DOF  are  computed  on  the  XDS-9300, 
which  also  does  the  magnitude  scaling  for  the  analog 
computer  (CI-500)  voltages.  In  addition  the  XDS-9300 
provides  the  required  control  logic,  timing  and  temporary 
storage  of  computed  variables. 

The  graphics  terminal  (AGT-10)  provides  real  time 
display  of  the  desired  output  and  control  variables,  as  well 
as  their  maximum  and  minimum  values.  Piloting  of  the  craft 
is  carried  out  at  the  graphics  terminal.  Additional 
permanent  record  outputs  of  navigation  x and  y are  obtained 
with  an  XY-pen  recorder,  and  plots  of  other  system  variables 
are  obtained  on  an  8 channel  Brush  recorder. 

A sample  of  the  x-y  output  for  the  3 DOF  simplified 
XR-3  craft  is  shown  on  Fig.  A- 3.  The  plot  starts  with  a 
right  turn  (5  = -0.5  radians)  at  t = 0,  using  a terminated 
ramp  input.  The  rudder  was  held  at  5 = -0.5  for  two  turns, 
then  reversed  to  5 = +0.5  for  a half  turn,  followed  by 
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5 =»  0.0,  then  again  3 = +0.5  for  a quarter  turn,  and  finally 
a right  turn  using  6 = -.29  radian  for  two  complete  turns. 

The  total  time  for  these  maneuvers  with  the  XR-3  was  approxi- 
mately 3^  minutes. 

Note  that  the  CI-5000  is  used  only  for  analog  inte- 
gration of  the  accelerations  to  obtain  the  desired  velocities 
for  the  digital  computation.  Displacement  is  computed  in 
the  digital  computer  using  the  trapezoidal  rule  for  integration 
as  done  on  the  Oceanic's  6 DOF  Loads  and  Motions  Program. 

Note  also  that  the  present  structure  and  programming 
use  only  existing  hardware  and  we  consider  them  a feasi- 
bility study.  If  results  show  that  we  have  been  successful, 
then  we  must  prepare  for  manual  control  studies,  and  this 
may  require  development  of  a more  realistic  "cockpit"  to 
simulate  conditions  in  the  control  room  of  the  3K-SES . 
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APPENDIX  B.  INSTALLATION  OF  THE  3K-SES , 5 DOF 
PROGRAM  IN  THE  IBM  360/67  COMPUTER 


In  the  fall  of  1978  the  3K-SES,  5 DOF  data  based 
program  called  DBSIM5D  was  delivered  to  the  Naval  Post- 
graduate School,  at  which  time  a commitment  was  made  by 
the  NPS  group  to  attempt  a hybrid  computer  simulation  based 
on  this  program. 

In  the  spring  of  1979,  DTNSRDC  shipped  the  card  deck 
and  the  output  data  for  a typical  3K-SES  run,  obtained  by 
DTNSRDC  on  a CDC  computer.  The  deck  was  to  be  used  in 
installing  the  program  in  the  IBM  360/67,  and  the  output 
data  provided  information  for  verification  and  validation. 
The  major  problem  in  the  process  of  installation  was  the 
conversion  from  CDC  Fortran  to  IBM  Fortran.  With  the  help 
of  Mr.  Roger  Hilleary  of  the  NPS  Computer  Facility,  the 
interactive  graphics  terminal  was  used  to  edit  the  program. 
After  some  50  man  hours  of  effort  the  first  successful 
batch  run  was  made  in  April  1979. 

Complete  success  in  the  conversion  did  not  occur  until 
a correction  was  added  to  the  newly  programmed  TIMER,  to 
allow  for  a stop  at  the  desired  stop  time  rather  than  a 
premature  stop  time  due  to  the  timer  program. 

To  validate  the  IBM  version  of  the  program  we  dupli- 
cated the  sample  run  for  the  3K-SES.  This  run,  as  provided 
by  DTNSRDC,  was  for  a left  turn,  using  an  11  degree  effector 
angle  and  a speed  of  60  knots.  The  sample  run  terminated 
after  14.5  seconds  of  real  (problem)  time.  The  11  degree 


l 

I 


r 


i- 


46 


effector  angle  was  programmed  as  a ramp  which  terminated 
after  1 second,  and  remained  constant  at  11  degrees  for 
am  additional  50  seconds. 

A comparison  of  the  output  of  the  IBM  360/67  version 
with  the  sample  run  data  from  the  CDC  computer  showed  that 
the  recorded  variables  were  in  complete  agreement  within 
the  accuracy  differences  of  the  two  computer  word  length. 

The  problems  of  conversion  were  simple  but  time 
consuming  and  may  be  listed  as  follows: 

1.  6 character  variable  name  limit  on  the  IBM  com- 
puter as  compared  with  a 7 character  limit  on  the 
CDC. 

2.  Single  statement  limit  on  the  IBM  computer  versus 
multiple  statement  on  a single  card  for  the  CDC. 

3.  Identification  of  temporary  files  used  for  the 
output  data. 

4.  Identification  of  the  timing  control  and  pro- 
gramming the  IBM  computer  timer  to  provide  the 
functions  required  for  operation  of  the  program. 

5.  Addition  of  the  proper  job  control  language  (JCL) 
cards  to  compile,  link  and  compute. 
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